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^-a  variational  technique,  an  approximate  dispersion  relation  is  found  for  arbitrary  density  profile,  and 
evaluated  in  closed  form  for  either  the  Bennett  or  square  profile.  Stability  properties  are  illustrated 
and  discussed  in  detail  for  a  square  profile,  including  the  influence  of  the  applied  magnetic  field 
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THEORY  OF  THE  RESISTIVE  HOSE  INSTABILITY  IN 
RELATIVISTIC  ELECTRON  BEAMS 


I .  Introduction 

There  is  a  large  and  growing  literature  on  equilibrium"  and  sta- 
2*6 

bility  of  intense  relativistic  electron  beams  propagating  in  initially 

T  8 

neutral  gas'  or  in  preionized  plasma.  The  most  serious  magnetohydro- 
iynamic  instability  of  a  beam  in  a  resistive  plasma  channel  appears  to 


be  the  resistive  hose  (m  =  l)  instability. ' 


115  Dacer  ieveicos  a 


Vlasov-Maxweil  theory  of  the  resistive  hose  instability  in  an  infinitely 
long  intense  relativistic  electron  beam,  propagating  parallel  to  a  uni- 

A  A 

form  magnetic  field  B  e  ,  with  axial  velocity  3-  ce  (Fi<z.  1).  A  dense 

o~z  c  ~z 

background  plasma  is  assumed  to  provide  space  charge  neutralization,  and 

.  2  2  2, 

the  beam  electron  motion  is  taken  to  be  paraxial  (p^  »  ?r  +  Pg~), 

2 

which  requires  that  v/y  «  1,  where  v  is  3udker's  parameter  and  y  me 

0  D 

is  the  characteristic  energy  of  the  beam  electrons.  Equilibrium  and 


stability  properties  are  calculated  for  rigid-rotor  equilibria,  which 


depend  on  energy  H  and  canonical  angular  momentum  P~  exclusively 
through  the  linear  combination  H  -  ui,  PQ  ,  where  uu,  is  a  constant,  and  for 

D  y  D 

the  related  cold  laminar  flow  equilibria,  where  ju^  can  be  a  function  of 

radial  position  r.  Some  of  the  results  presented  in  this  paper  appear 

o 

in  the  unpublished  literature,  but  are  included  as  required  for  com¬ 


pleteness  . 

Not*:  Manuscript  submitted  August  17,  1979. 
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Equilibrium  properties  are  examined  in  Sec.  II.  In  Sec.  IIIA,  an 

integro-differential  eigenvalue  equation  is  derived,  assuming  only  that 

k  R,  «  1,  juR,  «c,  and  that  there  is  complete  snace  charge  neutraliza- 
o  b 

tion,  where  k,  jj,  and  3^  are  the  axial  wavenumber,  frequency,  and  beam 
radius.  The  eigenvalue  equation  is  reduced  to  an  ordinary  differential 
equation  [Eq.  ( Uj ) j  in  either  the  cold  laminar  flow  limit,  or  in  the 
case  of  a  square  radial  density  profile.  3y  using  a  variational  tech¬ 
nique,  an  approximate  dispersion  relation  [Eq.  61)], 


is  obtained  in  Sec.  HID.  [A  slightly  modified  form  holds  when  there 

is  a  conducting  wall  at  R  .]  Here  d(r)  is  the  plasma  electrical  con- 

c  * 

ductivity,  ju,(t)  is  the  betatron  frequency,  uu  is  the  cyclotron  frequency 
P  c 

A 

associated,  with  the  field  3  e  ,  f  is  the  fractional  current  neutraliza- 

o 

A 

tion  of  the  equilibrium ,  and  A  (r)  e  and  J.  (r)  are  the  equilibrium 

O  DO 

vector  potential  and  beam  current  density,  respectively.  The  dispersion 
relation,  Eq.  (61),  constitutes  one  of  the  main  results  of  this  paper, 
and  can  be  used  to  investigate  the  resistive  hose  instability  for  a 
broad  range  of  system  parameters  and  arbitrary  beam  density  profile. 


In  the  special  case  of  a  Bennett  equilibrium  profile,  the  dispersion 
relation  can  be  further  simplified  to  ]Eq.  (76)] 


where  Tr  =  ^  is  a  characteristic  diffusion  time  [Eq.  (78)] 


for  the  perturbed  magnetic  field,  and  jj2  =  ju _(r  =  0)  is  the  maximum 

Pffl  fj 

value  of  ju,  .  This  is  a  generalization  (to  include  r.on-zero  f  and  ju  ) 
p  me 

of  a  result  obtained  by  LeeJ  within  the  framework  of  a  distributed- 
artificial-mass  model.  For  a  square  beam  profile,  the  dispersion  rela¬ 
tion  reduces  to 


-2 


U)a2  1-R,  2/R  2 
p  _ o  c _ 

1  O  0  5 

-m  1  -  uu  Td(l-R^  /Rc-) 


which  also  includes  the  effect  of  a  conducting  wall  at  R  . 

In  Sec.  IV,  the  stability  properties  are  illustrated  and  discussed 
at  length  for  a  square  beam  profile,  with  particular  emphasis  on  the 
influence  of  the  applied  field  (ju^),  the  conducting  guide  (through 
R^/F^),  and  the  fractional  current  neutralization  (f  ).  It  is  found 
that  the  applied  magnetic  field  and  the  finite  channel  radius  strongly 
reduce  the  growth  rate,  whereas  low  frequency  perturbations  are  con¬ 
siderably  destabilized  by  even  a  small  fractional  current  neutralization. 
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II.  Equilibrium 


We  consider  the  equilibrium  configuration  illustrated  in  Fig.  1, 
consisting  of  an  intense  relativistic  electron  beam  propagating  parallel 

A 

to  a  uniform  applied  magnetic  field  E_z  £  ,  through  a  background  plasma 

whose  conductivity  varies  with  radial  position.  Cylindrical  polar 

coordinates  (r,  9,  z)  are  used,  with  the  z-axis  along  the  axis  of 

symmetry.  Both  the  beam  and  plasma,  in  equilibrium,  are  assumed  to  be 

azimuthally  symmetric  (d/30  =  0),  infinitely  long,  and  uniform  axially 

(3/3z  =  0).  The  beam  is  characterized  by  density  profile  n.  (r) ,  cur- 

oc 

rent  density  J,  (r) ,  and  number  of  electrons  oer  unit  length  N,  and 

effective  radius  R_  ,  each  defined  by 
o 

00 

N.  =  2tt  f  drr  n.  (r)  »  n  3.  2  a.  (r  -  0) .  (l) 

0  Jr  DO  b  DO 

In  this  paper,  we  investigate  the  equilibrium  and  stability  proper¬ 
ties  of  steady  state  (3/3t  =  0)  beam  distributions  of  the  form 


r.  (H,  ?-,?)=  F  (K  -JU.P„)  5  ( P  -Pj, 
bo  0  z  b  a  zb 


where  the  total  energy. 


„  /  2  It  2  2, f  2 

H=(mc  +p  c  )  =  Y  me  , 


(2) 


(3) 


the  canonical  angular  momentum, 


P3  =  HPq  -  f  e3ozr/c)’ 


(U) 


and  the  axial  canonical  momentum. 


P  =  p  -  —  A  ( r 

Z  Z  C  0 


(5) 


are  the  three  single  particle  constants  of  the  motion.  In  the  above, 

£  =  (p  »  p0 i  p,)  is  the  mechanical  momentum,  and  -e  and  m  are  the  electron 

A 

charge  and  rest  mass.  A  (r)  =  A  (r)  e  is  the  vector  potential  associ- 

~  O  O  ~2 

ated  with  the  self-generated  azimuthal  magnetic  field  3~(r);  from 
Ampere's  law. 


r  —  A 
dr 


T  [Jbo(r) 
c  bo 


J  (r)l 

po 


f  )  J  (r), 
m  bo 


(6) 


where  J  (r)  is  the  axial  plasma  return  current.  For  simplicity,  we 
po  • 

consider  in  this  pater  only  the  case  where  J  is  a  ’uniform  constant 

po 

fraction  of  J-^tr),  as  indicated  in  Ea.  (6)_  We  note  that  the  distribu¬ 
tion  function  (2)  depends  on  H  and  ?a  exclusively  through  the  linear 
combination  H  -  D^Pg  ,  where  d.  is,  in  general,  a  constant.  (In  the  cold 
fluid  limit  discussed  below,  each  particle  is  localized  to  a  particular 

value  of  r,  and  it  is  Dossible  to  allow  id,  to  be  a  function  of  r,  thus 

b 

permitting  sheared  rotation.)  The  mean  azimuthal  motion  of  the  beam 
electrons  described  by  Eq..  (l)  corresponds  to  rigid  rotation  at  angular 
velocity  d^.  It  is  assumed  that  the  beam  is  completely  space  charge 
neutralized  by  the  background  plasma. 

We  assume  in  addition  that  the  motion  of  the  beam  electrons  is 
paraxial,  i.e. 


2  2  2 

p  +  p,  «  v  . 
rr  *  z 


(7) 
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Inequality  (7),  which  applies  to  the  random  component  of  electron 
motion  as  well  as  to  ordered  beam  rotation,  requires  that 


S,e'  (1  -  f  ) 


Y,  me 
o 


1" ,300  3.  Y-  amps 


« 


where  v  is  known  as  Sucker’s  parameter,  V,  =  3.„-  and  yvmc  are  the 

characteristic  beam  electron  axial  velocity  and  enersy,  and  I  ^  is  the 

—  met 

beam  current  less  the  plasma  return  current.  It  also  follows  from  (7) 
that  the  axial  self-generated  magnetic  field  'due  to  beam  rotation)  can 

be  neglected,  since  its  dynamical  influence  is  smaller  than  that  of  the 

o  o 

azimuthal  self- field  by  order  p.  ”/p_  .  Thus  the  equilibrium  magnetic 
field  is  oaken  to  be 


3  (r)  =  B  e  +  3  -,(r)  e  ^  , 

~o  oz  ~z  oa  ~  o  ’ 


where  3  ,  (r)  is  determined  bv  Eq.  (6). 

Using  (7),  it  follows  (after  some  straightforward  algebra; 
H  -  ju,  ?flcan  be  approximated  by 

C  0 


H  -  V9  *  Yb  mc"  +  2^  +  #0(r)’ 


where  _o  is  the  random  component  of  transverse  momentum. 


?/  5  P/  +  (?0-Vb^iib)2 


the  effective  potential  y  (r),  associated  with  both  magnetic  and 

0 

centrifugal  forces,  is 


2  r  2  2-1 

:l0(r)  =  f  Yb  mr  L  (juvju^v"  :  +  jj.  (r)J, 


1 


m  ~(r)  =  2ep,  A  (r)/(Yvmr2) 

fj  D  O  D 


(13; 


is  the  squared  betatron  frequency  (the  electron  oscillation  frequency 

in  the  absence  of  beam  rotation  and  of  an  axial  B-field) , 

ju  =  eB  /fflYtc  is  the  beam  electron  cyclotron  frequency,  and  the 
c  oz  b 

characteristic  axial  velocity  of  the  beam  is  defined  by 


V,  =  ,3,c  =  ?  /v.  a. 
0  b  0  0 


(it) 


In  (13)  we  have  set  A(r  =  0)  =  0,  without  loss  of  generality. 

Making  use  of  Eqs.  (2)  and  (10)  ,  the  beam  density  profile  can  be 
expressed  as 


.  X  „ 

n^Q(r)  =  2rr  C  dUE^-Unc  +  U), 

j  ill  (  -V 


(15) 


V  lr 
o 


2  . 


where  U  =  H  -  ju.  Pa  -  y.mc  is  an  effective  transverse  enerey  variable 
b  d  b 

in  the  rotating  frame  of  reference,  and  the  lower  limit  of  integra¬ 
tion,  )  (r),  is  the  minimum  possible  value  of  U  for  an  electron  at  r. 

2  2  2 

To  lowest  order  in  (p  +  p  )/?  ,  the  mean  axial  velocity  of  the  beam 

I*  j  z 

is  approximated  by  /v ,  defined  in  Eq.  (It),  and  is  independent  of  r. 


Similarly,  we  find  that  the  mean  azimuthal  velocity  is  approximately 

(16) 


Vr’  *  pr  VJ/ ^ 


ptbo  S  ^b1"  ’ 


corresponding  to  riaid  rotaticn  at  constant  angular  velocity  ju.  .  From 

0 

Eq.  (6),  and  the  constancy  of  Vv ,  we  obtain  the  equilibrium  azimuthal 
field 

Vr)  ■  -  -  f»)r“/0  drV  %o,r'  >• 


(1") 


From  Eqs.  (6),  (13)  and  (IT),  (r)  can  also  be  written 

2/  ,  l+TTe^  2  / 

3  - 2  (1  “  fm}  /  dr  r  nbo(  }< 

Y,  mr  •/o 

1  b 


(13) 


We  also  consider  the  case  of  cold  laminar  flow,  characterized  by 


i'll 


_ o 

dr 


(19) 


for  all  r.  Substituting  Eq.  (12)  into  Eq.  (19)  and  solving  for  w  , 

b 

ve  find 


ajv  =  uu  (r)  =  f  ju  ±  [juj  2  +  ju„2(r)]2  . 

D  C  C  p 


(20) 


The  equilibrium  described  by  Eqs.  (2)  and  (20)  is  one  in  which  each 
electron  is  localized  near  a  particular  value  of  r,  performing  a  nearly 
circular  helical  orbit.  Since  r  is  a  constant  of  the  motion  (to  lowest 
order  in  p  /y,m  r  uu  ),  uu  can  be  permitted  to  depend  on  r,  giving 

Xu  0  0 

sheared  cold  fluid  rotation,  rather  than  rigid  rotation.  The  conditions 
=  ju  (r)  and  uu„d  =  oj  (r)  correspond  respectively  to  fast  and  slow 
laminar  rotation. 

It  will  also  be  useful  to  rewrite  Eq.  (15)  in  the  form 
1  dnbo 


r  dr 


2tty.  m  Fp/mc  +  Ji  (r)l  — 
b  L  b  o  J  r 


1  a*o 


=  2rry  2m2F[y  me2  +  ty  (r)][ju  +  (r)-JJ,  ]  (r)j.  (21) 


As  a  first  example,  we  consider  a  loss  cone  distribution, 
VE'W  "  V2"  Yb«)  6(H^bP9-Y0BC2)6(P2-Pb), 


(22) 


3 


where  is  a  constant  satisfying  y„>Y-„- 


Eq.  (15),  we  find 
(r)  = 


00 


!n  0  <  r  <  R,  , 
= 

0  ,  R  <  r< 


Substituting  Eq.  (22)  into 


(23a) 


where  R^  is  defined  by 


2  2 

0  Y,  me  r  y  ~Y- 
_  2  c  o  o 

0  Yv 


(23b) 


and 


o  o 

V 

V 


(1  -  f  ) , 


(2U) 


so  that  uu  and  4  (r)/r'_  are  constants,  independent  of  r. 

H  O 

The  most  familiar  examples  of  beam  equilibria  with  diffuse  radial 
profiles  occur  when  the  beam  electrons  have  a  Gibbs  distribution  of 
"transverse  energy",  i.e. 

n. 


fbo(H’  ?9  ,?z }  =  2?T7^T  exp 


I  H^bpa-Yomc2  l 

\  H  '  )  5  <w  ■ 


(25) 


For  the  case  where  the  axial  magnetic  field  is  absent,  uu  =0,  or  where 

b 

its  effect  is  canceled  by  the  centrifugal  force  of  ordered  beam  rota¬ 
tion,  ju,  =  ju  ,  the  distribution  function  of  Eq.  (25)  generates  the 
DC’  w 

familiar  Bennett  profile.  Intense  relativistic  beam  equilibria  are 
discussed  in  many  references . 1 
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Ill .  Linearized  Vlasov-Maxwell  Analysis 


A.  Derivation  of  the  Integro-Differential  Eigenvalue  Sedation 
In  this  section,  ve  use  the  linearised  Vlasov -Maxwell  equa¬ 
tions  to  investigate  the  resistive  hose  (m  =  l)  instability  of  intense 
relativistic  beam  equilibria  described  by  2q.  (2).  We  adopt  a  normal 
mode  approach  in  which  all  perturbed  quantities  are  assumed  to  vary  with 
3  ,  z ,  and  t  as 

A 

i|f(x,  t)  =  'Jt(r)  exp  [  i (9  +  kz  -  Jut )]  ,  (26) 

with  Imuj>0.  In  practice,  it  is  often  more  convenient  to  use  T  and  z, 
rather  than  t  and  z ,  as  independent  variables ,  where 

T  =  t  -  z/V.  .  (27) 

D 

In  this  representation, 

I  (x ,t )  =  'If  ( r )  exp  [  i  (0  -  Gz  /V^  -  jut  )  ]  ,  ( 26  '  ) 

where  =  JU-kVfe  is  the  frequency  Doppler  shifted  to  the  beam  frame. 
Either  uu  or  D  may  be  regarded  as  an  externally  determined  quantity, 
depending  on  the  situation  envisioned.  If,  for  example,  a  perturbation 
is  initialized  by  a  deflecting  structure  at  z  =  0 ,  then  ju  is  the  real 
frequency  of  the  deflector,  and  Q ,  complex  in  general,  represents  the 
oscillation  and  growth  (or  damping)  of  the  wave  as  a  particular  beam 
segment  propagates  downstream.  If,  on  the  other  hand,  each  beam  segment 
is  taken  to  oscillate  at  a  fixed  real  frequency  fl ,  then  uu  represents 
the  growth  of  the  wave  as  one  moves  backward  from  the  head  of  the  beam. 

Q  occurs  through  the  dynamics  of  particular  beam  segments,  and  thus 


10 


scales  with  the  natural  oscillation  frequencies  of  the  bean  electrons, 
A-  ar.d  jj1  ,  while  a  occurs  through  the  magnetic  coupling  of  different 
beam  segments,  and  thus  scales  with  magnetic  decay  times. 

Ve  assume  that  the  plasma  is  coilisior.al  to  the  extent  that 
it  is  characterized  by  a  real,  scalar  conductivity  j(r).  Ve  further 
assume  that  the  perturbed  beam  space  charge  field  is  completely  neu¬ 
tralized  by  the  plasma,  which  requires  that 


!+tto  ( r ) »  a  , 


for  r £  R.  ,  the  beam  radius.  We  consider  only  wavelengths  long  and  fre¬ 
quencies  low  compared  to  quantities  that  characterize  the  beam  radius , 


A  ^  «  Cj 


(29b) 


therefore 


Q  Rj  «c. 


It  follows  that  the  transverse  electric  components  of  the  perturbed  fields, 

3  ,  E  ,  and  E- ,  can  be  neglected,  and  the  certurbation  can  be  represented 
z  r  3 

A  A 

in  terns  of  a  perturbed  axial  component  of  vector  ootential  A  e  : 

*■  *  *  'V 


BJr) 


(30a) 


3r(r)  =  p  A( r) , 


(3Cb) 


1  dA 
c  dt 


(30c) 


l: 
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Eqs.  (7)  and  (35),  the  eigenvalue  equation  (31)  for  resistive  hose  in¬ 
stability  can  be  expressed  as 


The  integro-dif ferential  eigenvalue  equation  (3 6)  constitutes  one  of  the 
main  results  of  this  paper  and  oar.  be  used  to  investigate  stability 
properties  for  a  broad  range  of  system  parameters.  We  emphasize  that 

Eq.  (36)  has  been  derived  with  no  a  priori  assumption  that  the  beam  is 

2 

cold,  or  'uniform,  or  that  the  perturbation  is  a  rigid  displacement ' 
However,  on  each  of  these  limiting  cases,  Eq.  (36)  can  be  simplified 
considerably,  as  discussed  in  the  remainder  of  this  section. 


3. 


Cola  Fluid  Description 


The  orbit  integral  I  and  corresponding  eigenvalue  equations 
are  obtained  in  this  section  for  the  cold  fluid  equilibrium. 
Substituting  the  identity 


oo  uu 

(r)  =  f  dr  '  v  '  Air  *)  f  dX  X  J,  (Xr  )  J  (Xr  '  ) 
Jo  Jo  i  j- 


into  Zq.  (36b),  we  have 

I  *  5T  /*  "  dr"  r"  Air")  [&  X  J.  (Xr"  )  f  °°id  7  exp  (  -ifl  T 
2n  Jo  Jo  1  ./  o 

2  it 

•  /*  d«S  J  (Xr)  exp  [i(97  -3'] 
yo  l 


(37) 


where  is  the  3essel  function  of  the  first  kine  of  order  l .  De¬ 
fining 


X  =  Xcos  a,  A  =  A  sin  a, 

X  y 


J.  (Xr  )  exp  [i(3/  -  3  )J  in  Eq.  (37)  can  be  expressed  in  terms  of  the 
Cartesian  orbit3  x  (r  )  =  r  (t  )  cos  3/  (t  )  and  y  (t  )  =  r  (T)sin  5*  (t  )  as 


J,  (Xr  )  exp  {  i(9'  -  9  ) } 


■/ 


2-rr 

da 

2TT 


—  exp  [i  (a  -  9-  '—)]  exp  (iX  x  +  iX  y  ) 

^  d.  X  J 


(38) 


In  the  cold  fluid  limiting  case ,  the  -unperturbed  transverse 
orbits  r  r  (%) ,  9*  (T)]  of  the  beam  electrons  are  near  circular,  i.e. 
r  *r.  It  follows  that  the  deviations  from  circiularity ,  of  order 

p>b2*V  (uu+  «i.  (39) 


J.4 


are  simple  harmonic,  and  the  unperturbed  orbits  are  given  by' 


X*  (7)  =  ~  ~~  C  3in(«i  +  ju+  7)  -  sin(  jS  +  ju~  7)  ] 

ju  -tu-  ‘b3 


+  r(uu, -uj-)  cos  (3-ku+7)  -  r(uu,-ju+)  cos  (9 -hju  7) 
b  b 


( UOa) 


y  (T  ) 


.-l- 

+  -  |  y,  m 

ju-ju  v  b 


: O S  (  35  +  JU  T  )  -  COS  \<£  +  JU  +  7)j 


+  r(ju,  -uu  )sin  (3  +  ju+ 7)  -  r(uu  -ju  +  )  sin  (9+ju  7! 

b  b 


I- 


( -Ob ) 


where  the  biharmonic  frequencies  ju  and  ju  ,  evaluated  at  r,  are  the  fast 
and  slow  laminar  rotation  frequencies  defined  in  Eq.  (20)  .  Substituting 
Eqs.  (Uo)  into  Eq.  (38)  and  carrying  out  the  integration  over  £,  we  have 

£  i£vXr')exp  ii(a'-9)J  I/0lex?CiM-?] 


{r(uu,  -ju”) 

O  r  , 

- + - - 

JU  -uu 


LX  cos  (3-hju+  7)  +  X  sin  (3tD+  T ) ]  (Ul) 
x  y 


r(uu  -ui.  ) 


[X  cos  (6-tui-  t)  +  X  sin  (9-hu-  t 

+  -  X  y 


JU  -0) 


“}• 


where  the  cold  fluid  approximation  (39)  has  been  invoked.  Simi-arly ,  the 

integration  over  cl  in  Eq.  (Ul)  gives 
2ir 

f  ^  J^(Xr)  exp  {  i(9'  -  9  ) } 


exp  ( iuu  T  )  + - exp  ( iuu 


X(r) 


(U2) 


i 
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where  the  effective  frequency  difference  x(t)  is  defined  by 


X(7  *f(  -U.  ~JU~f  +  (dJ+-U),  )2  +  (uu.  -dU~)  (0J+-U),  )  sin  (dJ^-dj")  7  J2  . 

0  DOC 


However,  the  cold  fluid  limiting  case  is  characterized  by  laminar  rota¬ 
tion^  at  either  oj~  or  jj+,  i.e. 


dj,  »  ju  ,  or  ju,  j)  ; 
o  b 


XCt  )  »  ju  -  ju 


Substituting  Sqs.  (1*2)  and  (1*5)  into  Eq.  (37),  and  carrying  out  the  in¬ 
tegration  over  time  and  r  /  ,  we  obtain 


_  + 

JU.  -dU  ,  JU  -dJ.  .  ”1 

o  x  _ b  - 

+  -  p.  +  +  ~  -  ’ 

m  JU  -UJ  i  JU  *\L)  ik 


which  is  independent  of  p^.  We  note  that  both  terms  in  (1*6)  are  essen¬ 
tial,  even  though  (UU)  indicates  that  one  of  the  terms  is  much  larger 
than  the  other;  the  larger  term  almost  cancels  out  when  inserted  in  (36a) 
The  integration  over  p^  in  Eq.  (36a)  is  now  trivial,  since  the 
term  in  brackets  is  independent  of  p  .  Using  Eq.  (21)  to  eliminate  7  in 

terns  of  J,  (r),  we  find  (after  some  algebra)  the  eigenvalue  equation 
bo 

for  the  cold  fluid  limit , 


±  i  A  rA(r)  ♦  maalli  A(r) 


A  (r) 


bo  (U7) 


Y,  mC  r  ~  r  \  7r  ^  ^  t  \ 

b  lO-jj  (r)j[Q-oj  (r) j 


+ ,  1 1  r  dr 


where  jf(r)  are  given  by  Eqs.  (20)  and  (l8),  in  terms  of  the  radial  pro¬ 
file  of  beam  density  n^Q(r)  or  current  density  J^Q(r).  The  treatment  of 


-o 


the  singularities  in  Eq.  (47)  is  determined  by  the  requirement  that  *1 
be  in  the  upper  half  plane.  The  boundary  conditions  on  Eq.  (k?)  are 

A.  (0)  =  A(Rc)  =  0,  <**8) 

where  may  be  30  if  the  beam  propagates  in  a  conducting  plasma  of  in¬ 
finite  extent,  or  is  finite  if  there  is  a  perfectly  conducting  guide  at 

r  =  R  .  Of  course,  the  case  of  finite  R  makes  sense  only  if  J,  (R  )=0. 
c  c  bo  c 

Equations  (U7)  and  (b8)  determine  the  hose  stability  properties  of  any 

beam  current  density  profile  J,  (r) ,  nrovided  that  the  beam  electrons 

bo 

are  in  a  cold  laminar  flow  equilibrium. 


In  Eq.  (36),  we  have  obtained  an  eigenvalue  equation  than  is 
very  generally  valid.  However  the  source  tern  in  Eq.  (36)  contains  an 

A 

integral  of  the  unknown  eigenfunction  A(r)  over  'unperturbed  orbits,  which 
makes  the  equation  rather  intractable  in  general.  This  difficulty  is 
fundamental,  reflecting  the  fact  that  individual  electron  orbits  span 
the  beam  cross-section,  communicating  information  about  the  perturbation 
from  one  value  of  r  to  another.  In  Sec.  III3,  we  avoided  the  problem 
by  considering  a  cold,  differentially  rotating  beam  equilibrium,  in 
which  each  individual  electron  is  confined  to  a  narrow  range  of  r;  Eq. 
(36)  then  reduced  to  the  ordinary  differential  equation  (1*7).  In  this 
section,  we  begin  to  explore  a  different  approach  in  which  the  orbit 
integral  is  evaluated  approximately  by  taking  advantage  of  simplifica- 

A 

tions  [periodic  orbits,  smoothly  varying  A(r)]  that  arise  naturally. 

In  this  section,  we  consider  only  the  case  of  a  beam  with  'uniform  den¬ 
sity  out  to  a  sharp  boundary  at  r  =  3^,  for  which  this  approach  is  most 
straightforward.  (More  general  density  profiles  will  be  considered  in 
subsequent  papers.)  We  do  not  assume  here  that  the  beam  is  cold;  never¬ 
theless,  for  the  uniform  density  beam,  we  arrive  at  the  same  eigenvalue 
equation  (U7).  For  generality,  we  continue  to  permit  the  presence  of  a 
guide  field  BQZ>  a  conducting  guide  at  r  =  (where  may  be  allowed 
to  go  to  ®),  and  an  equilibrium  return  current  density  f  . 

For  any  given  value  of  w  (which  plays  the  role  of  '&  in  the 
usual  three-dimensional  eigenvalue  problem) ,  there  is  in  general  an 
infinite  sequence  of  eigenvalues  ju  ,  corresponding  to  different  values 


A 

of  the  radial  quantum  number  [the  number  of  oscillations  in  A(r)].  Our 
interest  lies  in  the  lowest  eigenfunction,  which  corresponds  to  a  side¬ 
ways  displacement  of  the  whole  beam,  and  is  expected  to  be  the  most  un¬ 
stable  mode.  In  the  low  frequency  limit, 

UnaR^”  1  ju  j/c”«l,  (19) 

this  eigenfunction  becomes 

A 

A(r)  «  r,  (50) 

representing  a  rigid  displacement  of  the  unifora  beam  as  well  as  the 
self-generated  3q  field.  For  JJ  large  enough  that 

2  2 

IttgR^jj/c  Js  1, 

A 

magnetic  diffusion  begins  to  distort  the  eigenfunction  A(r),  even  if 

A 

the  beam  itself  were  to  displace  rigidly.  Nevertheless,  A(r)  follows 
(50)  closely  at  the  origin  and  over  much  of  the  beam,  bending  over  some¬ 
what  for  larger  r,  as  illustrated  in  some  examples  at  the  end  of  this 
section.  As  we  shall  see,  the  inequality  (49)  holds,  but  not  as  a 
strong  inequality,  over  the  whole  range  of  strongest  instability. 

In  view  of  this  discussion,  and  particularly  since  the  integral 

'  A 

in  Sq.  (36b)  is  merely  ar  average  over  A(r)  ,  and  is  not  expected  to  de- 

A 

pend  strongly  on  the  detailed  structure  of  A(r),  we  use  the  assumed 
form 

A(r)  =  A(r)r/r  (51) 

in  the  integrand  of  (36b).  Equation  (51)  is  correct  at  r  =  r,  as  well 


19 


as  having  the  right  form  at  the  origin,  and  becoming  exact  for  either 
small  r  or  small  uj  . 

For  a  uniform  beam,  ju?  and  ju~  are  independent  of  r,  and  it  is 
easy  to  show  that  the  unperturbed  orbits  are  given  exactly  by  the  simple 
harmonic  forms  of  Eq.  (lO).  Using  Eq.  (51),  Eq.  (36b)  reduces  to 


I  =  iA(r)r  ~  dr  exp  (-  i9-ift  T  )  f  (d«5/2rr)r  e~^ 

_co 

o  2tt 

=  iA(r)r  j_x  dt"  exp  (-i3+iH  7)  (i^/2rr)  +  i y  ), 


d~  exn  (-i3+if2  7) 


and  using  Eqs.  (10)  for  the  orbits  (x*,  y) ,  the  integrals  in  Eq.  (52)  are 
easily  performed.  The  result  is  identical  with  Eq.  (16).  Thus  Eq.  (IT) 
also  holds,  and  explicitly  inserting  the  'uniform  beam  current  density, 
we  arrive  at  the  eigenvalue  equation 


2  2 

ill  *,  N  ,  lnicr(r\ju  *  ,  N  ^ob  ^b 

- - —  r  Avr)  +  - ^ - A  (r)  =  - - - 

ar  r  ir  2  -,,0  +, 


5  ( r-R,  ) 


A(r),  (53) 


“ob  -  (W"/m  Yb'? 


is  the  beam  plasma  frequency;  we  note  that,  for  a  uniform  beam, 

Ju,2  =  S  Ju  2  3  2(1  -  f  ). 

3  pb  o  m 

The  appropriate  boundary  conditions  are  again  those  of  Eq.  (18).  «e  note 
that  if  Eq.  \19)  holds,  the  eigenfunction  A(r)  within  the  beam  is  indeed 
given  by  Eq.  (50) ,  as  discussed  previously. 

If  a(r)  is  a  constant  a  within  the  beam,  and  a  (possibly 
different)  constant  c.  outside  the  beam,  the  eigenfunctions  of  (53) 


are  Bessel  functions  of  order  one,  and  a  dispersion  relation  connecting 
x  and  n  is  easily  derived  in  closed  form,  from  the  boundary  conditions 
(1*3)  as  veil  as  the  matching  conditions  at  r  =  We  consider  the 

particular  case 


»  = 


j(  const)  ,  C<  r  <  R,  , 
o 

j,  <  <  c2/Utt  R,  2x  ,  R,  <  r  <  R 


The  condition  on  J,  in  Eq.  (5^0  insures  that  magnetic  diffusion  through 
the  weakly  conducting  region  F  <  r<R^  is  slow  compared  to  the  time 
scale  x  x ,  and  thus  that  terms  in  the  dispersion  relation  proportional 
to  3.  are  negligibly  small.  Strictly  speaking,  J  cannot  be  zero,  since 
our  assumption  of  electrostatic  neutrality  requires  that  0.  >  x  is 
easily  compatible  with  (5-*),  according  to  Eq.  '28).  Cases  in  which 
-rra  falls  below  ju  for  r  greater  than  some  radius  ?  can  easily  be  treated 

by  applying  appropriate  boundary  conditions  at  R  ,  as  discussed  by 

o 


-ee . 


For  a  beam  of  radius  R  ,  with  3(r)  given  by  Eq  ( 5L' 


A( r)  =  J. (Kr) 


,  0  ^  rsi, 


A(r)  =  f 


ci.  R  ^-r2 
b  c 


2  2 
R  -R, 
c  b 


V"  3^),  ^  ac. 


( 5o'o ) 


wnere 


,2  .  ,  2 

*  =  uni  <t<ji/c 


56) 


and  the  dispersion  relation  is 


(n  -  u»‘)(n  -  ju+)  =  x  .%.2f.(ub), 

pb  o  1 


(5”a) 


2 


t/'K 

ha 


D.  Variational  Approximation  to  the  Dispersion  Relation 

Except  in  special  cases  such  as  the  ’uniform  density  beam  with 
a  sharp  boundary,  considered  in  Sec.  Ill  C,  the  eigenvalues  and  eigen¬ 
functions  of  Eq.  (k?)  can  only  be  calculated  numerically.  Such  solutions 
will  be  presented  in  future  publications.  In  this  section,  we  develop 
a  variational  technique  of  approximating  the  eigenvalues  that  is  widely 


applicable . 


Beginning  with  Eq.  (1*7),  we  multiply  through  by  rA(r)  and  in¬ 


tegrate  over  r,  to  obtain 


IffTi  ju  r  2 

2  J 


A2 

drra(r)A  (r) 


R 

-  r  d 

'A  " 


irr 


i  A? 

lb  A2(r) 

^b  (Q-ou  )(D-uj  +  ) 


A(r)  A(r)  .  (58) 


)  .  ( ! 


If  we  guess  a  trial  form  A  (r)  for  the  eigenfunction,  and 


A"(r)  £  A(r)  +  6A(r ) , 


a 

where  A(r)  is  the  (unknown)  exact  solution,  let  uu  be  the  eigenvalue  cal- 
culated  (for  given  Q)  by  using  A  (r)  in  Eq.  (58).  We  also  let 


uu  =  uu  +  6uy 


where  ju  is  the  (unknown)  exact  eigenvalue.  We  find  then,  after  two  in¬ 
tegrations  by  parts,  that 


l*TTi5uu 


i6uu  f 

2  J  < 


A2 

dr  ra(r)  A  (r) 


-  f  R|rr«*(r)  [isaakl .  in.  fi. - i - i  fisa.  *  A  i  A  rl 

A  c2  c  %  (a-o<a-»+)  r  *  dr  r  ir  J 


5A(r),  (6i; 
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provided  that  the  trial  eigenfunction  A" ( r )  satisfies  the  required 
boundary  conditions  (U8).  Thus  the  first  variation 


and  the  calculated  eigenvalue  ju  is  accurate  to  second  order  in  the 


error  5 A. 


However,  the  differential  operator  in  Eq.  (17)  is  non- 


Hermitian,  because  Im  Q  >  0;  consequently  both  ju  and  A(r)  are  complex, 

and  the  condition  (62)  determines  a  saddle  point  in  the  complex  plane, 

rather  than  an  extremum,  as  in  the  usual  case  of  a  Hermitian  operator, 

where  the  eigenvalue  and  eigenfunction  are  real.  As  a  result,  the  error 

5uu  ,  although  small,  is  not  positive  definite,  as  would  be  true  for  a 

*  "fc 

Hermitian  operator,  i.e.  one  cannot  simply  seek  a  trial  function  A  that 
minimizes  juu.  Nevertheless,  Eq.  (6l)  is  valuable  if  one  can  find  a  trial 
eigenfunction  that  is  even  qualitatively  correct,  since  the  error  in  ju1' 
is  second  order.  An  obvious  choice  for  this  purpose  is  to  use  the  form 
of  A  (r)  corresponding  to  a  rigid  displacement  of  the  magnetic  field. 


A  (r)  «  dA  /dr, 
o 


provided  that  the  outer  boundary  condition  is  at  R  =  °°  .  [ If  R  is 

c  c 

At 

finite,  a  slightly  modified  choice  of  A  is  required  to  satisfy  the 

A 

boundary  condition  A(R  )  =  0,  Eq.  ( U8 ) ,  as  will  be  discussed.]  We  then 

find  that  Eq.  (58)  reduces  to  the  dispersion  relation 

.  r  00  /  dA  \  2  r  ®  dA  dJ.  /  uu„2(r)(l-f  )“ 1  \ 

—  f  drra(r)  (  )  =  f  irr  —^2.  ( l  +  - d - S -  |  ?  (61*) 

ojo  V  ar  dr  y  Cr^u"(r)][nV(r )}  J 


2k 


where  the  identity 


(C-uj  )(Q-0)+)  =  ft2-ou  il-jj  2 
c  p 


(65) 


has  been  used. 

We  consider  two  special  cases;  the  first  is  the  bean  with  a 

square  radial  density  profile,  also  discussed  in  Sec.  Ill  C.  For  this 

case  we  consider  also  the  presence  of  a  conducting  guide  at  R  .  The 

c 

natural  choice  for  a  trial  eigenfunction  is  then 


AL(r)  = 


V 


2  2 
R  -R 
=  b 


&-■). 


,  0^  r  2  R  . 

b 


(66) 


which  satisfies  (63)  for  r  <  R^,  but  net  for  r  >  R^,  and  also  satisfies 
(^8).  Equation  (53)  then  reduces  to  a  dispersion  relation  that  is 


slightly  different  from  (6u), 

2niw 


2a  2 
jj  3 

1  +  rb  ^b 


1-F^2/Rc2  (fUT)(Q-^  +  )  R  2c2 


*b 


I drr 


a(r ) 


(67) 
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R 

r  sl 

J  R.  r 


(H  2-r2)2 


Specializing  still  further  to  the  piece-wise  flat  conductivity  profile 
of  Eq.  (5M,  the  dispersion  relation  (67)  can  be  put  in  the  form 


with 


F2(uu)  =  - 


( 68b ) 


where  T  ,  essentially  a  decay  time  for  the  perturbed  current,  is  de¬ 
fined  for  this  case  by 


,  2.2 
Td  =  5  naR-o  /c  = 


2  2 

-  fi  *  \  Ml 


(69) 


and  K  is  the  quantity  defined  in  Eq.  (56).  Equation  (68a)  is  in 
exactly  the  same  form  as  the  exact  solution  of  ( 1?)  for  this  particular 
beam  and  conductivity  profile,  Eq.  (57),  except  that  an  approximate 
function  (ju  )  from  (63b)  replaces  the  exact  result  F,  (uu  )  from  (5Tc). 

To  check  on  the  accuracy  of  the  approximation  F^(dj ) ,  we  note  that 
F^U),  expanded  in  powers  of  juT^  up  to  second  order,  is 


F1(uu )  «  - 


1  +  (-  iUUT 


2 

)(1  -  H,2/R  2) 
b  c 


( 70 ) 


We  see  that  F^  and  F^  are  identical  up  to  a  numerically  small  correction 
2  2 

of  order  uu  ;  thus  the  dispersion  relation  of  Eq.  (68)  is  an  excellent 
approximation  for  jjt^  1.  This  is  as  might  be  expected,  since  the  model 
of  rigid  displacement  is  qualitatively  accurate  in  this  regime,  which  is 
expected  to  include  the  fastest  growing  modes. 

Alternatively,  Eqs.  (68a,  b)  can  be  written  in  the  forms 


JUT 


d 


-i 


1  -w 

(Q-uu  ~ )  (Q-ju+) 


(68c) 
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.  2  3  2 
=r2  2  ’ 

a  ^ 


(68d) 


where 


Cl2  =  Q(n-ai  ) 

C 


The  dependence  of  ju  on  Cl  is  given  explicitly  by  Eqs.  (68c,  d) .  We  note 
the  following  features:  (i)  For  w  real,  uu  is  pure  imaginary.  Instabil¬ 
ity  (Im  ju  >  0)  occurs  only  over  a  bounded  range  of  Cl,  which  is  a  subset 
of  the  range  ju  <  Cl  <  jj  +  .  (ii)  The  growth  rate  Im  ju  goes  to  infinity 
as  *1  -•  ju  from  above  or  Cl  -*  ju  from  below.  This  occurs  for  a  flat  beam 

profile,  because  all  beam  electrons  are  in  resonance  with  the  wave  at 

,  .  3 

these  frequencies.  For  rounded  beam  profiles  ,  the  growth  rate  is 

finite  for  any  value  of  Cl,  as  will  be  illustrated.  (iii)  Increasing  the 

2_2  2 

return  current  is  destabilizing,  since  uu  .  |3,  /uu_  =  2/(l-f  ).  Proximity 

pb  bp  m 

to  the  wall,  i.e.  increasing  R,  /R  ,  is  stabilizing.  The  effect  of  a  B 

DC  Z 

—2  2 

guide  field  is  merely  to  substitute  Cl  for  Q  ,  as  is  evident  from  Eqs. 

(1+7)  and  (65)  for  any  beam  profile.  Thus  the  growth  spectrum  Im  ju  ,  for 
Cl  on  the  real  line,  is  shifted  to  different  values  of  Cl  by  the  presence 
of  a  guide  field.  However,  we  shall  see  that  3^  is  stabilizing  when  Cl  is 
calculated  as  a  function  of  real  uu .  The  stability  properties  are  discussed 
further  in  Sec.  IV,  particularly  for  real  uu ,  and  with  emphasis  on  the 
effects  of  the  conducting  guide,  the  plasma  return  current,  and  the  3^  field. 
As  a  second  example,  we  consider  the  3ennett  profile, 
ntQ(r)  =  nbQ(0)(l  +  r2/!^2)"2  (72) 

which  is  the  equilibrium  profile  for  a  paraxial,  monoenergetic ,  non¬ 
rotating  beam  subject  to  isotropic  small  angle  scattering  by  a  uniform 


plasma  medium.  In  this  case,  the  dispersion  relation  (6U)  reduces  to 
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g  = 

0,  if 

2  „  pr2 
pm 

(79a) 

g  = 

0,  if 

.1^  <  0  (equivalent 

to  0  <  U  <  M  )  , 
c 

(79b) 

g  = 

-1,  if 

—2  2 

0  <  Q  <  uu_  ,  and  u  < 
pm 

0, 

(79c) 

g  = 

+1,  if 

—2  2 

0  <  Q  ,  and  0  < 

pm 

Q  (requires  ju  <  H). 

c 

(79d) 

For  the  special  case  with 


O 

II 

FJ 

Cm 

(80a) 

*c  =  °’ 

(80b) 

a(r)  =  ac  (1  +  r2/!^2)-", 

(31) 

Eqs.  (76)  and  (77)  reduce  to  the  dispersion  relation  previously  dis¬ 
cussed  in  detail  by  lee,"  who  arrived  at  this  result  by  using  an  en¬ 
tirely  different  set  of  plausibility  arguments — principally  a  model  in 
which  the  beam  is  sliced  into  overlapping  rigid  disks  with  an  artificial 


distribution  of  particle  mass. 

—2  2  —2 
The  function  G(Q  /<ju-  )  is  displayed  in  Fig.  2,  for  H  real, 

pm 

In  contrast  to  the  beam  with  square  radial  profile,  discussed  previously, 
the  growth  rate  Im  ju  is  bounded  (because  the  natural  oscillation  fre¬ 
quencies  vary  for  different  electrons,  and  there  is  no  single  fre¬ 


quency  Q  for  which  the  entire  beam  is  resonant  with  the  wave).  Another 

difference  is  that,  from  Eq.  (79)  >  Ju  is  complex  in  the  range 
—2  2 

0  <  71  <  .  Equation  (76a)  also  exhibits  explicitly  the  dependence 
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on  the  current  neutralization  factor  f 

m 

0  <  ftc  <  §  U0o  “is  unstable,  but  as  f  • 
pm  m 

— p  2 

-  00  <  ft-  ^  0 .65  JU~  becomes  '.instable. 

fj  m 


For  f  =  0 ,  only  the  range 
m 

,  the  entire  range 
According  to  Eq.  (70),  nega- 


— O  ,  P  _ 

tive  values  of  ii  up  to  -  j  JO  occur  for  ft  in  the  range  0  <  ft  <  ju. 


Growth  rates  Im  ju  also  go  to  infinity  as  f  -*  1.  The  destabilizing 

m 

effect  of  f  is  due  to  the  repulsive  interaction  between  the  beam  cur- 
m 

rent  and  the  return  current;  it  is  particularly  striking  and  significant 

that  instability  can  occur  even  for  ft  =  0,  when  f  >0. 

m 
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IV.  Effects  of  Guide  Field,  Conducting  Walls,  and  Current  Neutralization 
In  this  section,  we  present  an  analysis  of  the  hose  stability  prop¬ 
erties  of  the  beam  with  uniform  density  to  radius  R  ,  in  the  narrow  con¬ 
ducting  channel  of  Eq.  (5b).  The  dispersion  relation  is  plotted  in  a 
variety  of  ways,  to  illustrate  the  effects  of  the  applied  3  field,  the 

fractional  neutralization  of  the  eauilibriun  beam  current  (f  ),  and  the 

n 

conducting  guide  at  R  .  Some  of  the  results  in  this  section  are  similar 

9 

to  previous  work  reported  in  the  unpublished  literature,  but  are  pre¬ 
sented  here  for  completeness.  We  concentrate  on  the  region  jut  s;  ]_  0f 
parameter  space,  where  the  approximate  dispersion  relation  of  Eq.  (63), 
is  essentially  indistinguishable  from  the  exact  result,  Eq.  (57),  from 
which  results  are  plotted,  and  where  the  beam  displaces  essentially 
rigidly. 

Depending  upon  the  way  in  which  a  stability  problem  is  phrased 
(initial  value  problem  in  beam  or  laboratory  frame,  driven  oscillation) 
either  1  or  jj  may  be  regarded  as  the  independent  variable,  and  taken  to 
be  real;  for  general  ’understanding  of  convective  properties  of  the  wave, 
it  is  necessary  to  treat  both  variables  in  the  complex  plane.  However 
most  experiments  are  performed  by  "tickling"  the  beam,  at  real  fre¬ 
quency  ju  ,  at  a  particular  location  z.  In  Q  then  determines  the  growth 
of  the  wave,  on  a  particular  beam  segment,  as  the  beam  propagates  down¬ 
stream.  Our  discussion  here  will  thus  be  in  the  form  IUju),  where  ju  is 
real. 

We  note  again,  from  Eqs.  (57b),  (6b)  and  (65),  that  the  effect  of 

o  —2 

the  guide  field  is  to  replace  Cl~  with  Cl  a  Cl  {Cl  -  jj^)  in  the  dispersion 


3. 


relation.  Thus,  if  ou  is  regarded  as  a  function  of  real  G,  the  guide 
field  does  not  reduce  the  maximum  growth  rate  Im  ju  ,  hut  merely  shifts 
it  to  a  different  value  of  G.  On  the  other  hand,  for  a  given  real  value 
of  JJ ,  the  replacement  reduces  the  growth  rate  Im  Q,  as  will  be  amply 
illustrated,  and  thus  can  he  regarded  as  stabilizing  in  this  sense. 

We  consider  first  the  case  juT^~*  0,  i.e.  low  frequency  or  very 
high  resistivity.  The  dispersion  relation  then  reduces  to 


i<2-ju  g 


a.  A  2  +  fm 


CD  0 


0, 


(82) 


whence  a  necessary  and  sufficient  condition  for  instability  is 


f 

*m 


+ 


jj 


2_  2  ’ 
uu  3 
Db  pb 


(83) 


a  result  that  is  illustrated  in  Fig.  3.  It  is  indeed  remarkable  that 
instability  occurs  even  in  this  limit.  The  reason,  of  course,  is  that 
the  beam  tends  to  be  expelled  from  the  highly  conducting  central  region 
of  the  channel  by  the  return  current.  This  effect  is  particularly 
strong  when  the  highly  conducting  central  channel,  which  responds  to 
magnetic  perturbations,  is  narrow,  as  in  the  present  example.  The 
external  magnetic  field,  and  proximity  of  the  conducting  guide  at 
are  seen  to  be  stabilizing. 

To  further  illustrate  the  dependence  of  stability  properties 

on  the  fractional  current  neutralization  f  ,  nlots  of  (a)  the  growth 

m 

rate  .1,  and  (b)  the  real  frequency  G  ,  versus  irr  , ,  are  shown  in  Fig.  , 
i  r  d  D 
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for  ju  =  R,  /R  =0  ar.d  various  values  of  f  .  It  is  evident  that  as  the 
(.  o  c  m 

fractional  current  neutralization  increases  to  unity,  the  growth  rate 

increases  rapidly  (for  reasons  discussed  previously),  while  the  real 

oscillation  frequency  decreases  (because  the  magnetic  restoring  force  is 

being  canceled  out).  We  note  also  that  in  all  cases  the  maximum  growth 

rate  Q  occurs  for  jut  <  1,  which  is  the  reason  we  are  mainly  interested 

in  this  regime.  For  f  <<.1,  the  instability  is  driven  mainly  by  a  chase 

m 

lag  between  the  oscillation  of  the  beam  and  that  of  the  magnetic  re¬ 
storing  force,  and  the  growth  rate  peaks  for  JUT  ,  only  slightly  less  then 
unity,  but  for  f  >  0.*+,  instability  is  driven  mainly  by  the  simple 
repulsion  of  the  beam  current  and  the  equilibrium  current,  and  the 
growth  rate  is  remarkably  flat  for  jjt ^  ^  0.3,  actually  reaching  its 

maximum  at  jut  ,  =  0 . 

d 

In  Fig.  5,  ft.,  and  are  plotted  against  ju^  ,  for  f  =  0,  R  /R  =  O.h, 

and  several  small  values  of  jut  , .  We  note  that  increasing  the  guide  B 

d  oz 

field  (i.e.  ju  )  reduces  the  growth  rate  0.  considerably  (but  never  to 
c  i 

zero).  As  shown  in  Fig.  3,  as  well,  the  growth  rate  .1.  increases  with 

uu  T ^  in  this  range,  but  is  very  insensitive  to  jut 

The  influence  of  the  conducting  guide  at  illustrated  in  Fig.  6, 

where  H.  and  Q  are  plotted  against  uu  for  f  =  0,  jut  ,  =  0.5,  and 
l  r  c  m  d 

several  values  of  the  parameter  R^/R^.  As  shown,  the  growth  rate 
falls  rapidly  to  zero,  while  the  oscillation  frequency  increases,  as 
R^/R^  increases  to  ’unity.  The  reason  for  both  effects  is,  of  course, 
the  increasing  stiffness  of  the  magnetic  field  trapped  between  R^/R^. 
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In  Fig.  T,  H.  and  H  are  plotted  against  jj^  for  a  fully  current 

IT*  C 

neutralized  case,  f  =1,  with  R, /R  =  Q.U  and  jut  .  =  0.5.  Hie  growth 
21  1)  c  a 

rate,  due  to  repulsion  between  the  beam  and  the  opposite  plasma  current 

is  very  strong  for  ju  <  u  .  3,  ,  but  fails  off  when  ju  >  ju  ,3.  •  For  f  =  1, 

cpbb  cpbb  m 

u)_  -*  0  and  the  beam  is  essentially  confined  only  by  the  B  external 
3  oz 

field,  hence  the  linear  dependence  of  Q  on  -u  for  u  ^  ju  ,3.  ,  as  well 

c  roc  po  c 

as  the  strong  stabilizing  effect  of  increasing  ju  above  u  .3.  . 

c  pc  0 

We  summarize  this  section  by  noting  that  resistive  hose  instability 
can  be  driven  either  by  a  phase  lag  in  the  restoring  force  due  to  non¬ 
zero  jut  .  or  by  the  presence  of  a  plasma  return  current;  thus  increasing 
f  is  destabilizing.  Increasing  the  applied  field  ju^  reduces  the  growth 
rate,  but  never  to  zero.  Proximity  of  the  conducting  guide,  i.e. 

R^/R^  -  1,  also  exerts  a  strong  stabilizing  influence,  but  the  growth 

rate  remains  oositive  'until  R.  /R  =1.  These  general  features  of  the 

oc 

behavior  should  apply  to  a  beam  with  diffuse  profile  as  well  as  to  the 
•uniform  beam  discussed  in  this  section.  Additional  effects  of  great 
interest  occur  in  a  diffuse  beam,  due  to  mixing  of  electrons  with  diff¬ 
erent  resonant  frequencies.  These  effects  have  been  discussed  by  lee 
within  the  context  of  his  multiple-component  rigid  beam  model,'  and  will 


be  investigated  further  in  subsequent  papers,  using  the  formal  framework 


developed  in  this  paper. 


V 


Conclusions 


In  this  paper  we  have  formulated  a  Vlasov-Maxvell  theory  of  the  re¬ 
sistive  hose  instability  in  an  infinitely  long  relativisitc  electron 
beam  propagating  parallel  to  a  ’uniform  applied  axial  magnetic  field.  A 
dense,  resistive  background  plasma  provides  space  charge  neutralization. 
The  equilibrium  configuration  and  basic  assumptions  were  summarized  in 
Sec.  II.  A  formal  eigenvalue  equation  (36)  was  obtained  in  Sec.  Ill  A, 
and  was  reduced  to  an  ordinary  differential  equation  (47)  in  the  cold 
fluid  limit  in  Sec.  Ill  3,  and  in  an  approximate  treatment  of  the  specia 
case  of  a  square  beam  density  profile  in  Sec.  Ill  C.  Using  a  variationa 
technique,  an  approximate  dispersion  relation  was  obtained  in  Sec.  Ill  D 
for  arbitrary  beam  density  profile,  Zqs .  (64)  and  (6?  .  The  dispersion 
relation  was  evaluated  explicitly,  in  closed  form,  for  either  a  3ennett 
profile  or  a  square  profile,  including  the  effects  of  fractional  current 
neutralization,  the  applied  field,  and  (in  the  latter  case)  a  conducting 
outer  wall.  Stability  properties  were  illustrated  and  discussed  in 
Sec.  IV,  for  a  square  profile  over  a  broad  range  of  system  parameters. 

It  was  found  that  the  applied  magnetic  field  and  proximity  to  the  con¬ 
ducting  guide  reduce  the  growth  rate  of  the  resistive  hose  instability, 
whereas  even  a  small  amount  of  current  neutralization  is  strongly  desta¬ 
bilizing,  particularly  at  low  frequencies. 


The  authors  are  grateful  for  comments  by  and  discussions  with  Dr. 


William  Sharp. 
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stability  boundary  [Eq.  (83)]  in  parameter  space 
(w/wptA.  fm  -  Rb2/Rc2)  for  coTd  =  0 


